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■ To Pierre van Moerbeke on the occasion of his 60th birthday 

i Abstract. We give explicit formulas for Hadamard's coefficients in 

terms of the r-function of the KdV hierarchy. We show that some of 
the basic properties of these coefficients can be easily derived from these 
formulas. The first immediate corollary is the symmetry of Hadamard's 
coefficients about the diagonal. Another well known fact, which follows 
from this approach, is that on the diagonal Hadamard's coefficients de- 
termine the right-hand sides of the equations of the KdV hierarchy. 
^ ' The proof of the main result uses Sato theory and simple properties of 

\ Gegenbauer polynomials. 

o 
o 

in ', 1- Introduction and examples 



O 



Consider the one-dimensional Schrodinger (or Sturm-Liouville) operator 



Oh- d 2 

Its heat kernel H[x,y,t) is the fundamental solution of the heat equation 

d 



m i J / = a < L2 > 

^ . It is well known that H(t, x, y) has an asymptotic expansion of the form 

S3' 



H(x,y,t) ==- [l + y^ H n (x,y)t n ) as t -» + . (1.3) 

V J 

The differential equation (|1.2|) for H(x, y, t) implies the recursion-differential 
equations for the coefficients H n = H n (x,y): 

H = 1 (1.4) 
dH 

(x - y)-^ 1 + nH n = LH n ^i for n > 1. (1.5) 
ox 

This system is known to admit unique smooth solutions H n = H n (x,y) in 
some neighborhood of the diagonal x = y. The coefficients H n are named 
after J. Hadamard ^H], who constructed them for the first time. 

Computation of heat invariants of self-adjoint elliptic operators is a well 
known problem in spectral theory which has many applications, in particular 
to geometry and theoretical physics El EEH El El • The asymptotics 



Date: September 14, 2004. 



1 



2 



P. ILIEV 



of the one-dimensional Schrddinger operator are of particular interest due 
to their relations to the Korteweg-de Vries (KdV) hierarchy. More precisely, 
it is known that the restriction of the heat coefficients on the diagonal gives 
the right-hand sides of the KdV hierarchy, see |18[ I20| . 

In the present paper we show that there are simple formulas for the 
Hadamard's coefficients H n (x,y) in terms of the r-function of the KdV 
hierarchy (see the next section for a precise definition of the r-function) . 

A remarkable explicit formulas for the coefficients of the Taylor expan- 
sion of H n (x,y) around the diagonal x = y were previously constructed in 
[11 Theorem 1.3]. However, these formulas have a rather complicated com- 
binatorial structure and it is practically impossible to write a closed formula 
for the coefficients even for simple potentials u. One advantage of the for- 
mulas derived in this paper is that they give finite expressions for the heat 
coefficient if the r-function is known (e.g. the solitons, or the more general 
algebro-geometric solutions of KdV) . 

To see the importance of the KdV equations, let us compute the first 
few coefficients using the defining relations (|1,4|) - ()1.5|) , Anticipating the 
appearance of the r-function, let us write u{x) as 

u(x)=2 d2l °^ x) \ 
ox 2 

From (|1.5|) one can easily obtain simple formulas for Hi and H 2 

I™) (1 .6) 
x - y V t(x) T(y) ) 

and 

2 ' (% - y) 2 V V T ( x ) T (v) J 1 ' r(x)r(y) J ' 
For the third coefficient we have the following formula 

/ //// \ T "'( v \ 

(x-yfH 3 (x,y) = -6(x-y)H 2 (x,y) + 2l {) U ' 





t(x) r(y) 

2[ r"(x)r'(x) T"(y)r'(y) ^ | ^ r>{x)r»{y) r'(y)r"(xy 

Kt(x) r(y) r(y) t{x) / 



Notice that the integral cannot be computed explicitly, unless something 
remarkable happens. This is the place where the KdV equation comes in. 
Assume that u(x) depends on a additional parameter S3 and it satisfies the 
KdV equation 

4<9 3 w = u" + 6uu, (1.9) 
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where d% = d/ds^ stands for the partial derivative with respect to S3, and 
u 1 is the derivative with respect to x. Then one can easily see that 

' u 2 {m = _2 f r"'{x) _ r'"(y) \ + 2 (r"{x)r'{x) r"{y)r'(y) 



y 3 v T ( x ) T (y) J \ t2 ( x ) r2 (y) 

(1.10) 

which combined with Q1.8|) leads to a simple formula for H^(x,y). 

We'll extend these computations by showing that if u is a solution of the 
KdV hierarchy and r is the corresponding r-function, then there are simple 
explicit formulas for H n (x,y) in terms of r. 

The paper is organized as follows. In the next section we recall some basic 
facts about the KdV hierarchy and Sato theory, which are needed for the 
formulation and the proof of the main result. In Section 3, we prove a general 
formula for H n (x,y). It is interesting that the smoothness of the coefficient 
H n (x,y) on the diagonal is related to the Gegenbauer polynomials. As 
a corollary of the main theorem, we see the symmetry of the coefficients 
about the diagonal x = y as well as the connection between H n (x,x) and 
KdV equations. 

As another application of the explicit formula, we show in ^3] that the 
expansion is finite if and only if the potential u{x) is a rational solution of the 
KdV hierarchy decaying at infinity studied in [2] . Equivalently, one can 
characterize the corresponding operators as the rank one bispectral family in 
[Oj. For related results concerning the finiteness property of the heat kernel 
expansion on the integers and rational solutions of the Toda lattice hierarchy 
see ^2] ■ For solitons of the Toda lattice and purely discrete versions of the 
heat kernel see IE 



2. KORTEWEG-DE VRIES HIERARCHY AND SATO THEORY 

In this section we recall some basic facts about KdV hierarchy and Sato 
theory. For more details on this and the more general Kadomtsev-Petviashvili 
hierarchy we refer the reader to the papers |19l [7] or the more detailed ex- 
positions [51 OTj . 

Let 

d 2 

be a second order differential operator. The KdV hierarchy is defined by 
the Lax equations 

where j = 1, 3, 5, . . . is an odd positive integer and (L J / 2 ) + is the differential 
part of the pseudo-differential operator L^ 2 . The first equation (for j = 1) 
simply means that u(x, s\, S3, S5, . . . ) = u(x + si, S3, S5, . . . ), giving us the 
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convenience to occasionally identify x and s±. The next equation (for j = 3) 
is exactly the KdV equation (|1.9j) . Let us represent L in a dressing form 

L = Wd 2 W~ 1 , (2.2) 

where TV is a pseudo-differential operator of the form 

oo 

W = Y,^kd~\ Vo = l- (2.3) 

k=0 

The wave (Baker) function ^(x, s, z) and the adjoint wave function ty*(x, s, z) 
are defined as 



V{x,s,z) = TUexp \ xz + Y^ s 2i-i z2i 1 j 

(oo \ / oo v 

k=0 / V i=l 

and 



(2.4) 



<&*(x,s,z) = (W*) *exp I -xz - ^s 2i -i 



z 2i-l 



1=1 

oo 



(2.5) 



= ( fc ) exp I ~ xz ~ ^2 S2i ~ i * 

\k=0 / \ i=l / 

where W* is the formal adjoint to the pseudo-differential operator W. Using 
(|2.2jl one can easily see that 

L$>(x, s, z) = z 2 ^(x, s, z) and L^*(x, s, z) = z 2 ^*(x, s, z). (2.6) 

We shall also use the reduced wave function ^ and the reduced adjoint 
wave function ty* obtained from and respectively, by omitting the 
exponential factor, i.e. 

oo 

*(x,s,z) = Y J ^kZ~ k (2.7) 

k=0 

and 

oo 

^{x,s,z) = Y,r k z- k . (2.8) 

fc=0 

Equations (|2.6|) imply 

L^(x, s, z) + 2zd x ^(x, s, z) = and L^*(x, s, z) - 2zd x ^*(x, s, z) = 0. 

(2.9) 

Using equations (|2.2[) - (|2.5|) one can show that the wave and the adjoint wave 
function satisfy the following bilinear identities 



res- 



(z 2n ^«(x, s, z)V*{x, s, z)\ = 0, (2.10) 

for all nonnegative integers n and I, where ^ l \x, s, z) is the Zth derivative 
of \& with respect to x, and the residue is around z = oo. 
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The remarkable discovery of the Kyoto school was that the KdV hierarchy 
(|2.1j) could be described by a function t(x,s). The reduced wave and the 
reduced adjoint wave functions can be expressed in terms of t(x,s) by the 
following formulas 

r, \ t(x; s - [z^ 1 ]) - t(x;s + [z^ 1 ]) 

^f(x, s, z) = . — ^ — - and ^ (x, s, z) = . — ^ — -, (2.11) 
t(x,s) t(x,s) 

where [z] = (z, z 3 /3, z 5 /5, . . .). 

Finally, let us denote by W n (x,y) the coefficients of the function 1 
%(x,8,z)¥*(y,s,z), i.e. 

oo 

*(x,8,z)**(y,s,z) = Y J W n (x,y)z- n . (2.12) 

n=0 

Using 1)2.11)1 we can easily write an explicit formula for W n in terms of the 
r- function. If we denote by &k(s) the elementary Schur polynomials defined 
by 

oo / oo \ 

J2 &k{s)z k = exp (J2 "M-l**- 1 > (2-13) 

fc=0 \fc=l / 

then we have 

»fj. Jl ^» [6 ' |4| : M ; [6 :-' l *' l i (2.i4) 

T(x,s)T(y,s) 

where 

a = (8^83/3, d 2k ^/(2k - 1), . . . ). (2.i5) 

3. Explicit formulas for Hadamard's coefficients 

The main result of the paper is the following theorem. 

Theorem 3.1. The Hadamard's coefficients can be computed from the fol- 
lowing relation 

H n (x,y) = (-1)" E T ~ k{ n" ^ Wn - k ^ll (3-D 

where (a)k = oc(a + 1) . . . (a + k — 1) denotes the Pochhammer symbol, and 
W n (x,y) are defined by 1)2.14)1 . 

Proof. To prove that the Hadamard's coefficients are given by 1)3.1)1 we need 
to check that 1)1.5)1 holds and that H n (x,y) are smooth on the diagonal 
x = y. 

To see that 1)1.50 holds, let us denote 

/n( x,,,, ) = (-irE %;!: ) !ll fc ^- fc - 1 - (3-2) 

k=0 ' V y ' 

Then ()3.1)) can be rewritten as 

H n (x, y) = res 2 [f n (x, y, z)^(x, s, z)^*(y, s, z)] . (3.3) 
^This function is closely related to the Green function for L. 
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Using the last equation together with (|2.9|) one can easily see that 
[(x -y)d x +n) H n (x,y) - L(x,d x )H n _i(x,y) 

= res z ((x - y)d x f n {x, y, z) + nf n {x, y, z) - <9^/ n _i(x, y, z)) #(x, s, z)^*(y, s, z) 

+ {(% - y)f n {x,y,z) + 2zf n -. 1 (x,y,z) - 2d x f n - 1 (x,y,z))d x ^>{x,s,z)^*(y,s,z) 
A direct computation now shows that 

(x - y)d x f n (x,y,z) +nf n (x,y,z) - <9^/ n _i(x, y, z) = 
(x - y)f n {x,y,z) + 2zf n - 1 (x,y,z) - 2d x f n -i{x,y,z) = 0, 

which proves (|1,5|) . 

Next we need to show that H n (x, y) is well defined on the diagonal. Writ- 
ing H n (x,y) as 

2(-l)" ^ 2 fc (x - y) k {k + l) 2w _ 2fc _ 2 
ff »M=^pL (n-fc-l)! ^iH-l^y), (3-4) 

and applying L'Hopital's rule we see that we need to prove that for j = 
0, 1, . . . , 2n — 2 we have 



j., . . . , 

n-1 



g^(j)tf|ff^«(^)l^ = 0. (3.5) 
Using (f2~T2j) and 

<^§(x, s, z) = exp ^-xz - ^2 S2i~iz 2i ~^j {d x - z) j ^(x, s, z) 

we see that (|3.5j) is equivalent to the following identities 
fn-\ 



res. 







(3.6) 

Equation H3.6|) will follow from the bilinear identities H'2.1U[) if we can 
prove that the polynomial 

fc=0 v 7 v 1 

is an even/odd function when j is an even/odd number, respectively. It 
is a pleasant surprise to see that these polynomials are closely related to 
very well known classical orthogonal polynomials - the so called Gegenbauer 
polynomials. 

The Gegenbauer (or ultraspherical) polynomials are defined by 

^fsf^- 1 )^;^ >-!)', (3.8) 

see for example 3 , pages 302-303] . Notice that this definition can be used for 
arbitrary A. If A > —1/2 and A ^ these polynomials orthogonal on the the 



ON THE HEAT KERNEL AND THE KORTEWEG-DE VRIES HIERARCHY 7 

interval (—1, 1) with respect to (1 — x 2 ) x ~2 which, in particular, implies that 
C x (w) is an even/odd function when n is even/odd, respectively. However, 
we need these polynomials also for negative values of A. In this case, we can 
use the three term reccurence relation 

2(n + \)wC*(w) = (n + l)C x +l (w) + (n + 2A - l)C x _ x {w). (3.9) 

and Cq = 1 and C x = 2Xw to deduce that C x (w) is an even/odd polynomial 
when n is even/odd, respectively. 

Changing the summation index in Q3.7|) we can rewrite P n j(w) as 

fc=max(OJ-n+l) ' 

From the last equation and the defining relation (|3.8[) for the Gegenbauer 
polynomials one can see that 2 

P n j(w) = j\2 n - 1 (2n-2j-3)l\Cj~ 3 ~ s (w) for < j < n - 1, (3.11) 

Pn i(w) = ^77- C? J 2 M for n < j < 2n - 2, (3.12) 

JV ' (2j - 2n + 1)!! * v ; ~ > v ; 

which completes the proof. □ 

From l|2.14j) and (|3.1|) we obtain the following 

Corollary 3.2. The Hadamard's coefficients H n (x,y) are symmetric func- 
tions of x and y, i.e. we have 

H n (x,y) = H n (y,x). 

Finally, we show that the heat coefficients {H n (x, x)} determine the right- 
hand sides of KdV equations (|2.1|) . 

Corollary 3.3. We have 

2 n 

Hn(x, x) = ^W 2n (x, x) (3.13) 

(In — 1)!! 

and 

2n— 1 

[{L— )+,L] = 2d x W 2n {x,x). (3.14) 
Thus, the KdV hierarchy ()2.1j) is equivalent to the following equations 

(2n - 1)!! 

d 2n -iu= — 2w _ 1 " d x H n (x,x), for n = 1,2,... . (3.15) 



2 (-l)H = 1 and (2k - 1)!! = 1 • 3 • • • (2k - 1) for k > 1. 
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Proof. Using (j3.1j) and applying L'Hopital's rule 2n — 1 times we see that 
H n (x,x) = (-l) n 2 



(2n- 1)! 
n-l 

X 

k=0 

2 



^ L fc /2n-l\ (2n-fc-2)! 2w _ x _ fc 



(3.16) 

where P n ,2n-l{w) is the polynomial defined by (|3.1U|) for j = In — 1 

P„, 2 „_ lW . S 2— ( 2 "; 1 ) (3,7) 

Notice that this time we have 

Pn,2n-iH = (-l) n 2 n " 1 (2n - 2)!! (c^M + l) , (3-18) 

and using the same argument (the bilinear identity and the fact that C 2n _ 1 2 (if) 
is an odd polynomial) we obtain (|3.13j) from QH.lfijl . 
From ((221), (25), (ESI and (l2~TTjl it follows that 

L ^ = wd 2n-x w -l = &i{-d)T{x, S ) d 2n-l-i-j . gjCgM^ 

Combining this formula with (|2.14l) we see that the coefficient of d~ 1 in 

2n—l In — 1 

L — 2~~ is PF^Xjx). If we denote by (L - a - )_ the integral (Volterra) part of 

2n-l 

the pseudo-differential operator L 2 we obtain 

2n— 1 s , r , 2n-l . , , 2n-l , 2n-l . , , , 2n— 1 v , 

[(L— )+,L] = [(L—) + ,L] + = [L— - (L— )_,L] + = [L, (L~ )_] + 
= [9 2 +n(x),W 2n (x,x)a- 1 + 0(5- 2 )] + = 2^(^(3, a)), 
which gives (|3.14|) and completes the proof. □ 
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